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Performing shell model calculations for heavy nuclei is a long-standing problem in nu- 
clear physics. The shell model truncation in the configuration space is an unavoidable 
step. The Projected Shell Model (PSM) truncates the space under the guidance of the 
deformed mean-field solutions. This implies that the PSM uses a novel and efficient way 
to bridge the two conventional methods: the deformed mean-field approximations, which 
are widely applied to heavy nuclei but able to describe the physics only in the intrinsic 
frame, and the spherical shell model diagonalization method, which is most fundamental 
but feasible only for small systems. We discuss the basic philosophy in construction of 
the PSM (or generally PSM-like) approach. Several examples from the PSM calculations 
are presented. Astrophysical applications are emphasized. 

1. Introduction 

The nuclear shell model is the most fundamental way of describing many-nucleon 
systems fully quantum mechanically. However, using the conventional shell model 
based on a spherical basis to study deformed, heavy nuclei is a very difficult task 
because of large dimensionality and its related problems. Even with today's com- 
puter power, the best standard shell-model diagonalization can be done only in 
the full pf-shell space, as for instance demonstrated by Caurier et al. 1 , for which 
the dimension of the configuration space may already reach one billion. It seems 
impossible to apply this kind of shell model calculation to arbitrarily large systems. 

In performing calculations for a many-nucleon system, one should choose a 
proper shell-model basis to start with. While in principle, it does not matter how 
to prepare a model basis, it is important in practice to use the most efficient one. In 
this regard, we recognize the fact that except for a few lying in the vicinity of shell 
closures, most nuclei in the nuclear chart are deformed. This suggests strongly for 
shell model calculations to use a deformed basis to incorporate correlations in large 
systems efficiently. However, a deformed basis is associated with breaking of rota- 
tional symmetry and angular momentum in that basis is no longer a good quantum 
number. Although in some cases physics may be discussed in deformed bases at- 
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tached to the system (intrinsic frames), the broken rotational symmetry should, at 
least in principle, be restored by angular-momentum-projection, which corresponds 
to a transformation from the intrinsic frame back to the laboratory frame. Shell 
model diagonalization is then carried out in the projected basis defined in the lab- 
oratory frame. This is the scheme that the Projected Shell Model 2 is based on. It 
may represent a group of practical methods that allow one to perform shell model 
diagonalization calculations for deformed, heavy nuclei. 

In recent years, strong demand for a shell-model treatment for nuclei arises also 
from nuclear astrophysics. Since heavy elements are made in stellar evolution and 
explosions, nuclear physics and, in particular, nuclear structure far from stability, 
enters into the stellar modeling in a crucial way (see the recent review 3 by Apra- 
hamian, Langanke, and Wiescher). The nucleosynthesis and the correlated energy 
generation are not completely understood, and the origin of elements in the cosmos 
remains one of the unsolved physics puzzles. Nuclear shell-models can generate well- 
defined wave functions, allowing one to compute, without further approximations 
as often assumed in the mean-field methods, quantities such as transition probabil- 
ities, spectroscopic factors, and /?-decay and electron-capture rates for both ground 
state and excited states. These quantities provide valuable structure information for 
nuclear astrophysics. Indeed, in their example 4 , Langanke and Martinez-Pinedo has 
demonstrated that shell-model calculations could significantly modify the results of 
nuclear astrophysics. 

2. The projected shell model 

In performing shell-model-type calculations for large nuclear systems, the central 
issue has been how to truncate the shell-model space efficiently. It corresponds to 
a rearrangement of the configuration space to decouple the most important part 
from the rest of the space. There are many different schemes of truncating a shell- 
model space 5 ' 6, ' 8 ' 9 ' 10 . Although these schemes differ very much in details in the 
way of building the model bases and/or choosing effective interactions, they share 
the common character that the model space is first constructed by some physical 
guidance. The so obtained model space contain the most significant configurations, 
each of which can be a very complicated combination in terms of the spherical shell- 
model basis states. In this way, the basis dimension can be significantly reduced 
and the final diagonalization is carried out in a much smaller space, thus making a 
shell-model calculation for heavy nuclei possible. 

The present article reviews a particular truncation scheme implemented in the 
Projected Shell Model (PSM) 2 . One friendly feature of the PSM is its simplicity. It 
starts with the deformed Nilsson model combined by a BCS calculation to generate 
quasiparticle (qp) states. The Nilsson model 11 is one of the most successful single- 
particle models that has been carefully tested for different mass regions. In a PSM 
calculation, the shell model truncation is first achieved within the qp states with 
respect to the deformed Nilsson + BCS qp vacuum |0); then the broken rotational 
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symmetry (also broken particle-number and parity conservation, if necessary) is 
restored for these states by standard projection techniques 12 to form a shell model 
basis in the laboratory frame; finally a two-body Hamiltonian is diagonalized in 
this basis. The truncation obtained in this way is very efficient. For example, for 
excited rotational bands in a deformed, heavy nucleus (see examples below), quite 
satisfactory results can be obtained by a diagonalization within a dimension smaller 
than 100. In short, the PSM starts with deformed single particle states generated 
by (in principle any) mean-field solutions on one side, and performs a shell model 
diagonalization like a spherical shell model (but in a much smaller space) on the 
other. Such an approach lies conceptually between the two conventional methods: 
deformed mean-field model and spherical shell model, thus naturally bridges the 
two methods that have coexisted but rarely have any connections. 

The Nilsson + BCS calculation defines a set of deformed qp states (with 4 and 
4 being the creation operator for neutrons and protons, respectively) with respect 
to the qp vacuum |0). The PSM basis is constructed in the multi-qp states with the 
following forms 

. e-e: {|o>, 44lo), 44lo), 444410), 444410), 444410), . . .} 
. o-o: {4410), 444410), 4444|o), 44444410},...} 

• odd-^: {4|o),444lo)>444l°}>44444lo): • ■ •} 

• odd-7r: {4|o), 444l°)> 444l°)> 444441°). • • •} 

for even-even, odd-odd, odd-neutron, and odd-proton nuclei, respectively. A general 
term in above expressions may be written as \4> K ), with k denoting all the quantum 
numbers, such as the Nilsson quantum numbers, necessary to specify that term. 
The angular-momentum-projected multi-qp states are thus building blocks in the 
PSM wavefunction, which can be generally written as 

Wm°)= E ti° P MK\4>.)=Y.ti ap MK K \4>.)- (1) 

K,K<I K 

The index a labels states with same angular momentum and n the basis states. 
P MK is the angular-momentum-projection operator 12 and the coefficients f^' a are 
weights of the basis states. 

The weights j 1 ^" are determined by diagonalization of the Hamiltonian in the 
spaces spanned for various nuclear systems as listed above, which leads to the 
eigenvalue equation (for a given /) 

J2(H KK >-E a N KK ,)fZ,=0. (2) 

k' 

The Hamiltonian and the norm matrix elements in Eq. (2) are given as 

H KK . = (<f> K \HP^ Ki \4> K ,), N KK , = (<j> K \Pk K K> K , I'M- (3) 

Angular-momentum-projection on a multi-qp state \<p K ) with a sequence of / gen- 
erates a band. One may define the rotational energy of a band (band energy) using 
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the expectation values of the Hamiltonian with respect to the projected \<f> K ) 

E i = H^ = (^\HPk KK JM f4) 
K N KK (MPicKjM ' 

The central technical issue is how to compute the matrix elements, such as those 
in Eq. (3), in the projected states. This question applies generally to any models 
utilizing angular momentum projection. Following the pioneering work of Hara and 
Iwasaki 13 , a systematic derivation has been obtained for any one and two-body 
operators (of separable forces) with an arbitrary number of quasi-particles in the 
projected states 2,14 . In principle, the projected multi-qp basis recovers the full shell 
model space if all the quasi-particles in the valence space were considered in building 
the multi-qp states. However, the advantage of working with a deformed basis is 
that the selection of only a few quasi-particles near the Fermi surface is already 
sufficient to construct a good shell model space. The rest can simply be truncated 
out. 

In a usual approximation with independent quasiparticle motion (mean-field 
solutions), the energy for a multi-qp state is simply taken as the sum of those 
of single quasiparticles. This is the dominant term. The present theory modifies 
this quantity in the following two steps. First, the band energy defined in Eq. 
(4) introduces the correction brought by angular momentum projection and the 
two-body interactions, which accounts for the couplings between the rotating body 
and the quasiparticles in a quantum-mechanical way. Second, the corresponding 
rotational states are mixed in the subsequent procedure of solving the eigenvalue 
equation (2). The energies are thus further modified by the configuration mixing. 

If the deformed states have an axial symmetry, each of the basis states in (I), 
the projected \<fi K ), is a if-state. For example, an n-qp configuration gives rise to 
a multiplct of 2"~ 1 states, with the total K expressed by K = \K\ ± K<i ± • • • ± 
K n \, where Ki is for an individual neutron or proton. In this case, shell model 
diagonalization, i.e. solving the eigenvalue equation (2), is completely equivalent 
to a i-L-mixing. The amount of the mixing can be obtained from the resulting 
wavefunctions. 

The above discussion is independent of the choice of the two-body interac- 
tions in the Hamiltonian. In practical calculations, the PSM uses the pairing plus 
quadrupole-quadrupolc Hamiltonian (that has been known to be essential in nuclear 
structure calculations 15 ' 16 ) with inclusion of the quadrupole-pairing term 

H = H Q - l x ^QtQ, - G M P*P -GqY, PtK- ( 5 ) 

The strength of the quadrupole-quadrupole force \ is determined in such a way that 
it has a self-consistent relation with the quadrupole deformation e 2 • The monopolc- 
pairing force constants Gm are 



Gm — [Gi =F G*2 N a Z ] A 1 , 



(6) 
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with "— " for neutrons and "+" for protons, which reproduces the observed odd- 
even mass differences in a given mass region if G\ and G2 are properly chosen. Fi- 
nally, the strength Gq for quadrupole pairing is simply assumed to be proportional 
to Gm, with a proportionality constant in the range of 0.14 to 0.20, as commonly 
used in the PSM calculations 2 . 

3. Examples 

In this section, we show several examples from the recent PSM calculations. These 
include a light nucleus 48 Cr for which a full diagonalization of the pf-shell spherical 
shell model is available, a well-deformed nucleus 178 Hf in the heavy mass region 
with the well-known 31-year isomer, and a gateway nucleus 254 No to the superheavy 
"island of stability" . 

3.1. 48 Cr; A benchmark for models of pf-shell rotors 

48 Cr is a light nucleus for which an exact shell model diagonalization in the pf-shell 
model space and a mean-field analysis have been performed 17 . Though light, this 
nucleus exhibits remarkable high-spin phenomena usually observed in heavy nuclei: 
large deformation, typical rotational spectrum, and the backbending phenomenon 
in which the regular rotational band is disturbed by a sudden irregularity at a 
certain spin. In the PSM calculations 18 for 48 Cr, three major shells (N = 1, 2, 3) 
are used for both neutrons and protons. The shell model space is truncated at the 
deformation £2 — 0.25. 
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Fig. 1. Left panel: 7-ray energies _B 7 = E(I) — E(I — 2) in 48 Cr as functions of spin. Right panel: 
B(E2) values as functions of spin. These figures are redrawn from Figs. 1 and 2 of Ref. 18 . 



In the left panel of Fig. 1, the PSM result for 7-ray energy along the yrast 
band, together with that of the pf-shell model (pf-SM) reported in Ref. 17 , and that 
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of the Generator Coordinate Method (GCM) 18 , are compared with the experimen- 
tal data. One sees that the four curves are bunched together over the entire spin 
region, indicating an excellent agreement of the three calculations with each other, 
and with the data. The sudden drop in E 7 occurring around spin / = 10 and 12 
corresponds to the backbcnding in the yrast band of 48 Cr. In the right panel of 
Fig. 1, three theoretical results for B(E2) are compared with the data. All the three 
calculations use the same effective charges (0.5e for neutrons and 1.5e for protons). 
Again, one sees that the theoretical descriptions agree not only with each other 
but also with the data quite well. The B(E2) values decrease monotonously after 
spin 1 = 6 (where the first band crossing takes place in the PSM). This implies 
a monotonous decrease of the intrinsic Q-momcnt as a function of spin, reaching 
finally the spherical regime at higher spins. This implies also that the final results 
of a shell model calculation do not depend on choice of the basis (either spherical 
(£2 = 0) or deformed (e 2 = 0.25)); A spherical shell model can describe states with 
large deformation (low-spin states in Fig. 1) and a (projected) deformed shell model 
can describe states having spherical properties (high-spin states in Fig. 1). 

The results shown in Fig. 1 suggest that the PSM is an efficient shell model 
truncation scheme that reproduces the data with the similar quality as the large- 
scale pf-SM. To compare with the one-major-shell shell model 17 , the valence space 
in the PSM is much larger (three major shells) yet the shell-model basis is much 
smaller (< 100). An obvious advantage of the PSM is that because it works in 
a much smaller basis with the configurations having a clear physics meaning, the 
PSM is able to extract the undergoing physics easily. In another PSM calculation 19 
for the super-deformed band in 36 Ar, the same conclusion was drawn. The PSM 
reproduces the data in a much simpler way and yields a similar quality in the results 
as those of the large-scale shell model 20 . 

3.2. 178 Hf: A deformed heavy nucleus with a long-lived isomer 

Long-lived isomers - excited nuclear states with inhibited electromagnetic decay 
- may be considered to offer a form of energy storage 21,22 . The possibility to 
trigger the decay by the application of external electromagnetic radiation has at- 
tracted much interest and potentially could lead to the controlled release of nuclear 
energy 23 . The significantly high energy of the 178 Hf isomer leads to the expectation 
of a lower triggering threshold, on account of the higher level density. The high iso- 
mer energy (~ 2.5 MeV) and long half-life (31 years) in 178 Hf also leads to greater 
potential with regard to the utility of triggered energy release. In determining the 
favorable conditions for triggering 7-ray emission from the isomer, it is necessary 
to have information on the structure of possible gateway states, as well as possible 
paths of electromagnetic transitions to and from these states. Experimentally, only 
a few states close to the 178 Hf isomer have been observed 24 ' 25 . 

In the PSM calculation 26 for 178 Hf, the model basis is built with the deformation 
parameters e 2 = 0.251 and £4 = 0.056. The valence space includes three major shells 
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(N = 4, 5, and 6 for neutrons and N = 3, 4, and 5 for protons). Fig. 2 shows the 
calculated energy levels in 178 Hf, compared with the known data 24 . Satisfactory 
agreement is achieved for most of the states, except that for the bandhead of the 
first 8~ band and the 14" band, the theoretical values are too low. The leading 
structure of each band can be read from the wave functions. We found that 
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These states, together with many other states (not shown in Fig. 2) obtained from 
the same diagonalization, form a complete spectrum including the high- if isomeric 
states and candidate gateway states. 
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Fig. 2. Comparison of calculated energy levels in 178 Hf with known data 24 . This figure is adopted 
from Ref. 26 . 



This example has demonstrated that the PSM can be an appropriate nuclear 
structure theory for studying high-spin isomers and associated excitations, including 
potential gateway states for triggering the isomer decay. Since the diagonalization 
mixes the if-states, the resulting wavefunctions contain valuable information on 
JsT-mixing. One is thus able to use the wavefunctions to calculate the inter-band 
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transitions and to analyze the degree of if -mixing. 

3.3. 254 No: A gateway nucleus to the superheavy "island of 
stability" 

One key question in our understanding of the chemical building blocks of the Uni- 
verse is what is the heaviest element that can exist. Intimately linked to this question 
is the problem in nuclear physics of the shell stabilization of superheavy nuclei. In 
a simplistic model of the nucleus, the repulsive force between more than 100 posi- 
tively charged protons is sufficient to overcome the attractive strong nuclear force 
and induce the nucleus to fission. However, protons and neutrons fall into shell 
structures which, in analogy to the enhanced stability of noble gases, give extra 
binding to nuclei in a closed shell configuration. The largest binding comes from a 
coincidence of closed shells for both protons and neutrons, as in 16 (Z = N = 8), 
40 Ca (Z = N = 20), 132 Sn (Z = 50, N = 82) and 208 Pb (Z = 82, N = 126). 

What is the next closed shell for protons, i.e. where is a region of long-lived, or 
even stable, superheavy elements beyond the actinides? Modern theoretical calcu- 
lations disagree on the size and position of the island of stability. For the possible 
sources that lead to the disagreement, one may question the single-particle levels 
in those calculations. One accessible way to gain information on these states with 
today's experimental situation is to study the high-if isomer states in Z « 100 
nuclei. These isomer states are highly significant as their structure is very sensitive 
to single-particle levels. Recently, Xu et al. 27 have suggested that the existence of 
the isomer states can result in increased survival probabilities of superheavy nuclei. 

The PSM is employed to study the structure of high- if isomer states in 254 No 
and the obtained results are compared with data from the recent experiment by 
Herzberg et al. 28 . In the PSM calculation, three major shells are activated for both 
neutrons and protons (N = 5, 6, 7 for neutrons and 4, 5, 6 for protons). The 
deformed basis for the 254 No calculation is constructed with Si — 0.25. For the 
pairing interaction strengths in the Hamiltonian, we take G\ = 21.24, Gi — 13.13, 
and G Q = 0.14. 

We first compare the theoretical results with the known 254 No data 28 for the 
yrast band. As can be seen in the left panel of Fig. 3, the theory has well reproduced 
the data. The current data ends at I = 20 where the rotational alignment process is 
about to start. After that spin, the theory predicts a weak up-bending in moment 
of inertia. In the right panel of Fig. 3, band energies (as defined in Eq. (4)) of 
several important configurations are plotted to understand why this up-bending 
occurs. While the 0-qp band lies low in energy till I = 20, the \ow-K 2-qp band 
{K" = l+,i/[743]f~ ©i/[734]§~) crosses the 0-qp band at I — 22 and becomes the 
lowest one for the next spin states. At I = 30, the 2-qp band is crossed by the 
low-.fr 4-qp band {K* = 0+, f[743]|~ © *4 734 ]f~ © tt[633]| + © tt[624]§ + ), and the 
latter becomes the lowest one for the high-spin region. The predicted up-bending 
effect centering at I = 26 in the left panel is thus attributed to the two successive 
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Fig. 3. Left panel: Comparison of the predicted yrast band in 254 No with known data 28 . Right 
panel: Selected band energies before band-mixing, which represent the most important configura- 
tions. 

band crossings with presence of strong band interactions. 

Among the high- if 2-qp states, the one with K n = 8~ coupled from two neu- 
trons (t/[624] | + © z/[734]| ) is found to be the lowest in energy. The next lowest 
8~ state is from two protons (7r[514]| ©7r[624]| + ). These two 8~ bands are found 
well-separated from the rest 2-qp bands. The predicted bandhead energies for the 
two 8~ bands are 1.266 and 1.516 MeV, respectively, very close to the experimental 
8" isomer at 1.293 MeV 28 . 

The combination of these two 2-qp 8~ states can give a 4-qp state with the 
highest possible K. This K 71 = 16 + state thus has the configuration (z/[624] | + © 
^[734]§~ ©tt[514]|~ ©tt[624]| + ). This is the lowest 4-qp band found in the calcu- 
lation, which is again well-separated (in energy and in K) from other 4-qp states. 
Although the bandhead state (I = 16) already lies in a dense region, its unique 
K quantum number hinders its decay. The predicted bandhead energy for the 16 + 
band is 2.750 MeV. Very likely, this is the structure for the experimental 184-/zs 
isomer 28 . 

4. Application in nuclear astrophysics 

Nuclear structure input is important, and sometimes crucial, for nuclear astro- 
physics. For example, an isomeric state can communicate with its ground state 
through thermal excitations in hot astrophysical environments. This could alter 
the elemental abundances produced in nucleosynthesis. We are just beginning 22 to 
look at the impact that isomers have on various nucleosynthesis processes such as 
the rapid proton capture process thought to take place on the accretion disks of bi- 
nary neutron stars. Shell model calculations for understanding the isomer structures 
and the decay rates are very much desired for this study. 
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4.1. Nuclear shape isomers in the rp-process path 

It has been suggested that in x-ray binaries, nuclei are synthesized via the rapid 
proton capture process (rp process) 29 , a sequence of proton captures and /3-decays 
responsible for the burning of hydrogen into heavier elements. The rp process pro- 
ceeds through an exotic mass region with N « Z, where the nuclei exhibit unusual 
structure properties. One of those is the coexistence of two or more stable shapes in 
a nucleus at comparable excitation energies in nuclei with A w 70 - 80. The nuclear 
shapes include, among others, prolate and oblate deformations. In an even-even 
nucleus, the lowest state with a prolate or an oblate shape has quantum numbers 
j7r _ q+ An excited + state may decay to the ground + state via an electric 
monopole (E0) transition. For lower excitation energies, the E0 transition is usu- 
ally slow, and thus the excited + state becomes a "shape isomer" . 
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Fig. 4. Energy surfaces for various spin states in 68 Se and 72 Kr as a function of deformation 
variable £2- This figure is adopted from Ref. 30 . 



Fig. 4 shows calculated total energies 30 as a function of the deformation variable 
£2 for different spin states in 68 Se and 72 Kr. The configuration space and the inter- 
action strengths in the Hamiltonian can be found in the previous calculations for 
the same mass region 31 . Under these calculation conditions, it is found that in both 
nuclei, the ground state takes an oblate shape with £2 ~ —0.25. As spin increases, 
the oblate minimum moves gradually to £2 ~ —0.3. Another local minimum with 
a prolate shape (e 2 ~ 0.4) is found to be 1.1 MeV ( 68 Se) and 0.7 MeV ( 72 Kr) high 
in excitation. Bouchez et al. 32 observed the 671 keV shape-isomer in 72 Kr with 
half-life r — 38 ± 3 ns. The one in 68 Se is our prediction 30 , awaiting experimental 
confirmation. Similar isomer states have also been calculated by Kaneko, Hasegawa, 
and Mizusaki 33 . 

Since the ground states of 73 Rb and 69 Br are bound with respect to the isomers 
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in 72 Kr and 68 Se, proton capture on these isomers may lead to additional strong 
feeding of the 73 Rb(p, 7) 74 Sr and 69 Br(p, 7) 70 Kr reactions. However, whether these 
branches have any significance depends on the associated nuclear structure param- 
eters, such as 

• how strong is the feeding of the isomer states? 

• what is the lifetime of the isomer with respect to 7-decay and also to /3-decay? 

• what are the lifetimes of the proton unbound 69 Br and 73 Rb isotopes in compar- 
ison to the proton capture on these states? 

The lifetime of the isomeric states must be sufficiently long to allow proton 
capture to take place. No information is presently available about the lifetime of 
the 68 Se isomer. Based on Hauser-Feshbach estimates 29 the lifetime against proton 
capture is in the range of as 100 ns to 10 /is, depending on the density in the 
environment. Considering the uncertainties in the present estimates a fair fraction 
may be leaking out of the 68 Se, 72 Kr equilibrium abundances towards higher masses. 




Fig. 5. Mass fractions in the x-ray burst model with two extreme cases. This figure is adopted 
from Ref. 30 . 



While it is likely that equilibrium is ensued between all these configurations 
within the presently given experimental limits a considerable flow towards higher 
masses through the isomer branch cannot be excluded. Figure 5 shows the compar- 
ison between the two extreme possibilities for the reaction sequence calculated in 
the framework of a multi-mass-zone x-ray burst model 34 . The left panel shows the 
mass fractions of 64 Ge, 68 Se, and 72 Kr as a function of time neglecting any possible 
isomer contribution to the flow. The right panel shows the results from the same 
model assuming full reaction flow through the isomeric states in 68 Se and 72 Kr 
rather than through the respective ground states. The main differences in 68 Se and 
72 Kr mass fractions are due to rapid initial depletion in the early cooling phase of 
the burst. This initial decline is compensated subsequently by decay feeding from 
the long lived Ge abundance. The results of our model calculations are based 
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on upper and lower limit assumptions about the role of the shape isomer states. 
Improved calculations would require better nuclear structure data to identify more 
stringent limits on the associated reaction and decay rate predictions. 

4.2. Weak interaction rates 

The knowledge on weak interaction processes is another important ingredient for 
resolving astrophysical problems. The decisive role played by the Gamow- Teller 
(GT) transitions has since long been recognized. It has been suggested that the 
nuclear shell model is the most preferable method for GT transition calculations. 
This has recently been emphasized by Langanke and Martinez-Pinedo 4 . 

For a theoretical model employed in GT transition calculations, it is generally 
required that the model can reproduce a wide range of structure properties of 
relevant nuclei. It has been shown that the conventional shell-model diagonalization 
method is indeed capable of performing such calculations. For example, Wildenthal 
and Brown 35 ' 36 obtained nuclear wave functions in the full sd-shell model space, 
which were successfully applied to calculation of GT rates in the sd shell nuclei 37 . 
Langanke and Martinez-Pinedo 38 made the shell-model GT rates available also 
for the pf shell nuclei. Still, these sophisticated calculations are tractable only for 
nuclei up to the mass-60 region, and cannot be applied to heavier nuclei which play 
important roles in the nuclear processes in massive stars. 

In GT transition calculations for heavy nuclei, the PSM wavefunctions can be 
used. Here, we wish to list a few attractive features in this approach, which may be 
relevant for future astrophysical applications. 

• Because of the way the PSM constructs its basis, the dimension of the model 
space is small (usually in the range of 10 2 — 10 4 ). With this size of basis, a 
state-by-state evaluation of GT transition rates is computationally feasible. This 
feature is important because in stellar environments with finite temperatures, 
the usual situation is that the thermal population of excited states in a parent 
nucleus sets up connections to many states in a daughter by the GT operator. 
Our current knowledge on GT transitions from excited nuclear states is however 
very poor, and in many cases, it must rely on theoretical calculations. 

• The PSM is a multi-shell shell model. This feature is desired for the processes in 
which forbidden transitions are dominated by the collective response of nuclei. 
The corresponding calculations must involve nuclear transitions between different 
harmonic-oscillator shells, and thus require a shell model working in a multi-shell 
model space. 

• The existence of isomeric states in nuclei could alter significantly the elemental 
abundances produced in nucleosynthesis. There are cases in which an isomer of 
sufficiently long lifetime can change the paths of reactions taking place and lead 
to a different set of elemental abundances. We have seen in above discussions 
that the PSM is indeed capable of describing the detailed structure of isomeric 
states. 



9, 2008 0:42 WSPC/INSTRUCTION FILE nucl'ysun 



Shell model for heavy nuclei and its application in nuclear astrophysics 13 

Very recently, we have developed a method 39 for calculation of GT transition 
rates based on the PSM. We expect that this method can generally be applied to the 
fields where weak interaction processes take place in nuclear systems. In particular, 
one may find interesting applications to cases where a laboratory measurement for 
certain weak interaction rates is difficult and where the conventional shell model 
calculations are not feasible. Potential applications in nuclear astrophysics are cal- 
culations of /3-decay rates for the r-process and the rp-process nucleosynthesis, and 
electron-capture rates for the core collapse supernova modelling. 

5. Summary 

It has been desired that one enjoys advanced shell model diagonalization methods 
while pushing the calculations to heavy nuclear systems. However, most of the 
updated spherical shell model calculations are severely limited in the fp-shcll nuclei 
or nuclei in the vicinity of shell closures. We have shown that the Projected Shell 
Model may be an efficient truncation scheme for a shell model solution. While it 
takes the advantages from the mean-filed results, the PSM also contains the shell 
model character that it performs configuration mixing in a truncated basis. We have 
shown through the calculation for 48 Cr that the PSM results are comparable with 
those obtained by the large-scale shell model calculations based on the spherical 
basis. We have also presented examples for heavy and superheavy nuclei. 

The use of a schematic interaction makes the theory extremely powerful for 
practical applications. It is rather amazing that such simple schematic forces can 
describe a large number of different types of nuclei consistently and accurately. It 
may imply that they indeed represent the most important parts of the effective 
nuclcon correlations in nuclei, as discussed in Refs. 15,16 . 

An important shell model application is to provide the structure information and 
decay rates for nuclear astrophysics. There have been some shell model applications 
along this line (see, for example, Ref. ); however, it still has a very long way to go. 
A large, unexplored region is the region away from the /3-stability, in which nuclear 
astrophysics is much interested. There are many questions for every shell model, 
including the PSM, to answer; the immediate ones are those concerning the choice 
of proper effective interactions and single particle states for the new regions. 
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